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ABSTRACT 

The measurement of the integrated optical polarization of weakly gravitationally 
lensed galaxies can provide considerable constraints on lens models. The method out- 
lined depends on fact that the orientation of the direction of optical polarization is not 
affected by weak gravitational lensing. The angle between the semi-major axis of the 
imaged galaxy and the direction of integrated optical polarization thus informs one of 
the distortion produced by the gravitational lensing. Although the method depends 
on the polarimetric measurement of faint galaxies, large telescopes and improved tech- 
niques should make such measurements possible in the near future. 
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1 INTRODUCTION 

Gravitational lensing of distant galaxies by foreground galaxy clusters provides a powerful method of determining the distri- 
bution of matter on a large scale. Its strength lies in the fact that it is sensitive to all gravitating matter, and that it gives 
the matter distribution directly, without invoking other physics or other assumptions. In this way it is different from other 
methods based on such me chanisms as the Sunyaev-Zeldovich effect which depend on the on the hydrodynamical state of the 
cluster (( Wu fc Fang 199tj ) and references therein). 

There are of course weaknesses of the method which arise from the limited data available, and the difficulty with which 
these data are obtained. With the improvement of technique we can expect these problems to be less of an obstacle in the 
future. 



Considerab le effort has gone into i nv erting images pro duced by gravitational lensing by clusters ((Kochanek 1990), (Lupine 
Kaiser 1997), (Bonnet et al. 1993), (seitz et al. 1996)) to find the properties of the gravitational field and mass of the 



lensing cluster. One distinguishes two regimes - that of strong lensing, where the mass density is near the critical mass density 
for the lensing system, which produces arcs, extreme distortions, and multiple images of galaxies, and weak lensing where 
only minor distortions are observed. 

In this paper we concentrate entirely on the weak lensing regime which is appropriate for the outer part of galaxy cluster. 
In the weak lensing regime one is looking for distortions of background galaxies which are at high red-shift and very faint. 
Redshifts are difficult to determine for most imaged galaxies, and the determination of their shape subject to considerable 
uncertainty. The effect of weak lensing on an elliptical source, as might be provided by an spiral galaxy inclined to the line of 
sight, is to change its apparent size, orientation and ellipticity. In principle measurement of this change can yield information 
about the gravitational field of the lens, and hence the matter distribution within the cluster. In practice this is difficult because 
one lacks information. Even if one assumes that both the distance of the lens and that of the source can be determined from 
redshift analyses, this still leaves an uncertainty in the intrinsic size, ellipticity and orientation of the source galaxy. These 
uncertainties have led researchers to treat the problem as a statistical one, where the backgr ound galaxies are tre ated as a 
statistical ensemble without any attempt to use information about their intrinsic properties. (( Bonnet et al. 1993| ), (Kaiser et 
al. 19951), (ISeitz et al. 19961)) 



If in some way one's ignorance of the source parameters could be reduced the lens parameters could be better determined. 
Measurement of the apparent brightness and redshift are obvious ways to improve inferences, but are difficult to achieve given 
the faintness of the images. Another interesting property that has yet to be exploited is the optical polarization of the source 
galaxy. Measurement of galactic polarization of high redshift galaxies would of course be difficult, given the fai ntness of the 
objects, although this has bee n carried out for a number of radio galaxies in both radio ( Gabuzda et al. 1992 ) and optical 
bands (Tadhunter et al. 1992). Potentially it could yield valuable information. (Indeed, in the case of radio observations 
these measurements have been carried out with considerable accuracy with a resolution of milliarcseconds.) The reason is 
the following. Light polarization is not affected, at least if the lensing object has small angular momentum (in relativistic 
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units) (schcinder et al. 1992) by gravitational lensing in the sense that the direction ol polarization will not be altered. For 
a typical spiral galaxy containing free electrons and dust in the galactic disc the integrated optical polarization arising from 
s cattered starlight c ould be as high as 2 percent and one would expect it to be in a direction perpendicular to the disc 
(Bianchi et al. 1996). Even if other mechanisms such as dichroic extinction are at play, the reasonable assumption of large 
scale axial symmetry requires that the integrated polarization should be either in the direction of the major axis of the ellipse 
or perpendicular to it. Thus the integrated polarization could in principle be used to fix the orientation of the source galaxy, 
modulo 90 degrees. 

The degree of polarization will obviously depend on the angle of inclination of the source galaxy to the line of sight and 
the number of scattering particles in the galactic disc. I f we assume Thomson o r Rayleigh scattering and axial symmetry in 
the distribution of sources and scatterers, one can show ( Simmons & Audit 1995 ) that the degree of polarization for optically 
thin galaxies is proportional to sin 2 i and the optical depth along the galactic radius, where i is the angle of inclination of the 
axis of symmetry of the galaxy to the line of sight. 

One signature of a lensed galaxy would be an observed deviation of the direction of polarization from the perpendicular 
to the semi-major axis of the image ellipse. In principle it is possible for the direction of integrated polarizati on to be in th e 
direction of the sem i- major axis o f the ellipse, owing to othe r mechanisms as discussed by several authors (( Draper 1995 ), 
(scarrott et al. 199C), (Wood 1997) and ( Wood and Jonesl997)). Small scale dev iations from rotational symmetry produced by 
dumpiness in the scattering distribution, as discussed by ( |Witt fc Gordon 1996| ), or inhomogeneities in the galactic magnetic 
field in the case of dichroism, should not greatly affect the direction of integrated polarized flux although the degree of 
integrated polarization could be reduced. 

In the ideal situation the degree of polarization should yield information about the inclination of the source galaxy to 
the line of sight. It would in fact, however, be difficult to accurately measure the degree of polarization, or to unambiguously 
infer from it the inclination, and hence the ellipticity on the sky of the spiral galaxy. (High inclination gives rise to high 
polarization.) 

In this paper we investigate how measurement of galactic light polarization can be used to sharpen and improve the 
determination of the lens parameters. In section (Q) we give a rough estimate of the degree of polarization produced in a spiral 
galaxy due to dust and electron scattering. Our main purpose in doing this is to ob tain a ball-park value fo r the degree of 
polarization. More accurate modelling has been carrie d out for the optically thin c ase (|Simmons fc Audit 1998 ), and by Monte 
Carlo treatment of radiative transfer ((Wood 1997), (Wood and Joncsl997) and (Bianchi et al. 1996)) for the optically thick 
case, but with the view to obtaining a resolved polarization map rather than the integrated polarization. In section we 
discuss how information on the inclination of the source galaxy together with measurements of the ellipticity and orientation 
of the images can be used to infer the parameters of the lens, and in section (^]) we compare the results obtained when this 
polarimetric information is available with the case when it is not, in order to assess the usefulness of the method. The last 
section contains a short discussion of the feasibility of carrying out such a polarimetric measurement and our conclusions. 



2 GALACTIC POLARIZATION 



Measurement of optic al polarization of spiral galaxies has been carried out for a few nearby galaxies ((scarrott et al. 1990) 



(Ncininger et al.l99C)). Various mechanisms could produce this optical polarization, the most obvious being scattering by 
dust and electrons. Most recent studies however have largely been concerned with the detailed structure of the galaxies and 
the possible existence of appreciable magnetic fields. Accordingly they have aimed at providing a resolved map of polarization 
in optical and infrared, which could give some information not only about the distribution of scattering grains and electrons, 
but also about the presence of magnetic fields and of dust lanes near the galactic nucleus. Dust grains, assumed oriented 
by the magnetic field of the galaxy, preferentially absorb light a long their longer axis , and hence could produce polarization 



perpendicular to this axis of the grain (dichroism). Such a model (Wood and Jonesl997) has been partly succesful in explaining 
the observed pattern for several high inclination spiral galaxies where the extinction would be most pronounced. In such cases, 
near to t he nucleus the polarizatio n is parallel to the galactic plane. However, that is not necessarily the case for other spiral 



galaxies (Scarrott & Draper 1996) 



In this paper we are interested only in the polarized flux integrated over the whole galaxy, which should in principle give 
an indication of the orientation of the galaxy. If the spiral galaxy has an axis of symmetry, so that its isophotes are ellipses, 
we should, on the grounds of symmetry, expect the linear polarization flux to be either perpendicular the major axis of the 
ellipse, or possibly, if such effects as dichroism were important, along the minor axis. Although the small scale deviations 
from axial symmetry have been observed for spiral galaxies at high inclinations, we should expect most of these effects to be 
largely cancelled out in the integrated flux, which we would expect to be either perpendicular to the major axis of the ellipse 
or along the major axis. The apparent rotation brought about by weak lensing would not usually be more than a few tens 
of degrees, and thus the determination of the direction of integrated polarization 'modulo' 90 degrees would be sufficient to 
provide pretty unambiguous information about the orientation of the source galaxy. 

In the case of spiral galaxies, one would expect optical polarization to arise as a result of Thomson, Rayleigh and 
dust scattering of starlight primarily in the galactic disc of the spiral galaxy, with minor contributions from the halo. For 
Thomson and Rayleigh scattering, light scattered through 90 degrees will be entirely polarized in a direction perpendicular 
to the scattering plane. (For other scattering mechanisms we would expect the polarization to be either in the plane or 
perpendicular.) Thus if Thomson or Rayleigh scattering were dominant, a galaxy inclined to the line of sight should display 
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Observer 



Figure 1. Simple scattering geometry. The galaxy is represented by a annulus, with inner radius Ri and outer radius R2, with the 
source at the centre. The angle between the line of sight and the axis of symmetry (X) is denoted i. The ray scattered by element at a 
is indicated by dot-dashed line. 



polarization, and one would expect in the case o f a rotationally symmetric disc that the direction of polarization to be along 
the minor axis of the ellipse (Bianchi et al. 1996), although in exceptional circumstances it could be along the major axis. 



2.1 A simple model for the integrated polarization 

One can show that under fairly general assumptions the degree of polarization produced by Thomson or Rayleigh sc attering, 
in the single scattering regime depends on sin 2 i and the total number of scattering particles ( [gimmons fc Audit 1996 ). In this 
section we shall give a very simplified derivation of this result to obtain an order of magnitude for the galactic polarization. 

Consider a completely flat galactic disc with the light source at the centre (see Fig. [jj ). Take the surface density of 
electrons in the disc to be Af e and the luminosity of the galaxy to be L. The axis of symmetry of the galaxy, which we have 
taken to be the z-axis, is at inclination i to the line of sight. We take the x-axis to be in the plane of symmetry. 

The flux arriving at the scattering element at a is given by L/4iva 2 . Thus the energy scattered per steradian per unit 
time into the line of sight by this element is simply 

La 



dF = 



, „AfJa)adOda—- (1 + cos x) 



where 9 is the angle between a and the x-axis, \ is the scattering angle and a the total scattering cross section. We obtain 
similar expressions for the polarized flux. In terms of the Stokes parameters referred to the scattering plane we have 



dF Q = -r—^Afe(a)ad6da-^- sin 2 % 

47TGT l07T 



and 



dF v = 0. 

To obtain the total scattered energy per steradian per unit time we integrate over the disc to obtain 

Z" 7 , Afe(a)a-^—(l + cos 2 x)ddda. 
47ta^ lD7t 

Using the properties of the Stokes parameters under rotations, total polarized energy Fq expressed in the observers frame 
(with the polarimeter aligned with the axis of symmetry) is given by 



„-A/'c(a)a sin 2 y cos 2d> dOda 

iira z 167T 
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and 

Fu = I / ^ a o A/" e M11-7- sin 2 x sin 20 rfflrfa , 
J J 47rar lD7r 

where is the angle between the scattering plane and the plane defined by the line of sight and the axis of symmetry of the 
galaxy. The scattering angle x an d the angle <j> may be expressed in terms of 9 and i 

cos x = — cos 9 sin i 

and 

cos 9 cos i 



sin x 

Substitution of these into the integral expressions yields after a bit of reduction 

La 3n . 2 . f J\f e (a) 

bn = sin 1 I da 

y 4tt 8 J a 

and 

Fu = 0, 

signifying that the polarization lies along the axis of symmetry (because Fq is always greater or equal to zero). We shall 
ignore the unpolarized scattered power as this will be swamped by the direct light from the galaxy. The degree of polarization 
is simply Fq j direct wnere ^direct ^ s tne direct power per steradian, L/4ir. Thus the degree of polarization is 

3tt a . 2 . f Af e (a) , 

p = sin 1 I da . 

8 J a 

The degree of polarization depends on sin 2 i. There are obvious problems in evaluating this integral for most density distribu- 
tions because of the possible singularity at a — which arises from the two dimensional treatment of the problem. To avoid 
such difficulties we take the surface density to be zero within a certain radius of the galactic centre. Again since we are really 
only interested in orders of magnitude here, we take the surface number density to be constant within an annulus with inner 
radius Ri and outer radius i?2- We then obtain for the degree of polarization 

37rA/"eO- . 2 • , R2 

p = sin 1 In ■ 

^8 ili 

If we denote the optical depth through the galactic disc by r and assume a geometric thickness for the galactic disc of A 
then the degree of polarization may be written 

3T7T . 2 . A R 2 

p = sin 1— —In-—-. 

8 R2 — -Ri Ri 

Typically we can expect the ratio of the thickness of the disc to its radius to be 0.1. so that the degree of polarization is 
approximately p = O.lr sin 2 i where we take the value of InRijRx to be of the order of 1. Thus with an optical depth in the 
range 0.1 to 1 we obtain a maximum polarization of one to ten percent. The actual polarization depends of course also on the 
inclination of the galaxy through the sin 2 i factor. 

This simple mod el seems to give the right order of magnitude for the integrated polarization when compared to more 



sophisticated models ( Bianchi et al. 199€| ) 



3 THE LENS EQUATION AND WEAK LENSING 

A simplified picture of the lensing situation is represented in figure ^| which serves to define the quantities entering into our 
equations. The lens equation relates the position, 77 of the source galaxy (in the source plane which is a distance Ds from 
the observer) to the position, £, of the image galaxy in the lens plane, which is a distance to the observer. (Distance here 
means the angular size distance). Elementary geometrical arguments yield 

n = TT-Z - ^ds"' (!) 
^d 
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Source plane 



Lens plane 

Figure 2. Simplified picture of the lens situation. 



Observer 



where a' is the vector angular deflection produced by the lens. If we use angular coordinates instead this equation may be 
written 



y = x - a, 



(2) 



where a is defined to be a'D a s /D s- The deflection angle is determined by the matter distribution in the lens plane. With 
the usual assumptions (see Schneider et al. 1992) we may write 



where 



a = V* 



^(x) = — / d x' k (x')\n\x — x'\ 

TV 



(3) 



(4) 



and k' = E/E cr ^, and £ cr ;t = c 2 Ds /A-kGD^ D^ s , where E is the projected mass density. It follows that 

V 2 * = 2k' . 

The amplification matrix, A, is then defined as the inverse of the jacobian matrix (dy l / dx J ), and so yields 



1 - *n *12 
*12 1 - *22 



(5) 



It follows from a linear expansion o f equation H tha t in the weak regime, an elliptical source (i.e. elliptical isophotes) will 
be transformed into an elliptical image (Kochanek 199C). The relation between the two depends on the properties of the lens, 
and is given below by equations (^-^) in terms of the of the characteristics of the source and image ellipses, and the lens. 
Following Kochanek, we characterise the source ellipse by the triplet S = (X s , A\ s ,a s ), the image by X = (Xi, A\i,cti) (a 
gives the orientation of the ellipse and (A + AA) 1//2 (resp. (A — AA) 1//2 ) its major (resp. minor) axis length) and the lense by 
£ = (k, 7, 9). Effectively these provide 



A\ a S a = (k - 7 2 )S l AA l (6) 

A s = A,( K 2 + 7 2 ) +2AX,K 1 C t (7) 

A S ±AA S C* S = (A<± AA l C* l )(ft± 7 ) 2 (8) 

AA S C* S = 2k7A 1 + AA 1 C* 1 (k 2 + 7 2 ) (9) 

where we have defined: C' x — cos2(aa; — 6*) and S x = sin 2(a x — 9). There are 9 unknown parameters and three independent 
equations. Thus in principle any one of S, X,or C can be determined from a knowledge of the other two. In particular the 
source parameters, S, can be determined if one knows the image, X, and lens parameters, C. 

In what follows we shall also use the parameter e describing the source and image ellipticity defined in terms of A and 
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AX as 

'A- AA 



e= 1 



A + AA 



It is interesting to note, before making any attempt to solve these equations, that they have a number of scaling 
laws which should make their solution easier. Thus for a given lens, if (Ai, AAi, on, X s , AX S , a s ) is a solution, then so is 
(xXi, xAXi, <Xi,xX a ,xAX a , ot s ) where a; is a constant. Moreover, for a given image, if (A s , AA S , a s , k, 7, 6) is a solution, then so 
also is (y 2 X s , y 2 AX s , ot s , yK, yj, 6). This means that one can arbitrarily take A^ and k to be 1. Dividing the equations 
by n 2 Xi (x — 1/AAi and y — 1/k) one obtains: 

Al s Ss = (l-r 2 )SiAk (10) 

Is = (1 + T 2 ) + 2AUTC l (11) 

l s ±Al s C's = (1± AZ 1 C l )(l±r) 2 (12) 

AhC a = 2r + A/iCi(l + r 2 ) (13) 

where Ah = AXi/Xi, l s = X s /(K 2 Xi), Al s = AX 3 /(K 2 Xi), and F = j/k. 



4 INFERRING THE LENS PARAMETERS 

The determination of the lens parameters in the weak lensing regime is usually treated as a statistical problem. In the absence 
of additional information about the source galaxies, in particular about the orientation of their axis of symmetry (if they have 
one), most analyses to date have attempted to solve the problem of inferring the lens parameters by matching the distribution 
of source parameters obtained from the data by the lens model, with an assumed distribution of source galaxy parameters e.g. 
a uniform distribution of orientations of source galaxies. One can obtain the best fit model param eters by standard routines 



such as x fitting. This sort of approach has been adopted by , for example, (Bonnet et al. 1993). Although very promising, 
the measurements are not good enough to provide accurate determination of the lens parameters. The techniques developed 
to reconstruct the cluster potential from the le nsed galaxies are n o w extremly sop his ticated and the y treat with great care 



the observational and statistical uncertainties ((Kaiser et al. 1995), (Schneiderl996), (Kochanek 199C)). We do not intend in 



this paper to go into such details. Our point is just to show that polarimetric measurements put strong constraints on the 
cluster potential reconstuction. 

We have already argued that information about the source galaxies, and in particular about their orientation on the 
sky, which could be obtained from measurement of their optical polarization, would lead to more accurate determinations. 
Polarization measurements would undoubtedly be difficult for such faint galaxies. However, even the measurement of the po- 
larization for a small subset of the galaxy sample can considerably improve the accuracy on the lens parameters determination. 
In the following section we illustrate this numerically for simple lens models. 

To assess what can be obtained via polarization measurements, we have carried out a simple simulation. First we decide 
on a lens model, which fixes the values of k(x, y), y(x,y) and 9(x,y) in the lens plane. We then randomly generate a sample 
of source galaxies from a given distribution with the following properties. We take the distribution of the symmetry axis to 
the line of sight to be isotropic, which implies a uniform distribution in orientation on the sky, and uniform distribution in 
ellipticity (if we assume the galactic discs to be circular). Since we are going to use the system which contains only 

"reduced" variables, we do not need to know the brightness of the galaxies at this stage. These galaxies are then lensed by 
a given lens model to give a distribution of images. We then attempt to reestablish the lens parameters from the sample of 
images. We shall see below that the effectiveness of the retrieval depends crucially on the assumptions about the model, and 
whether or not we use potential information about the orientation of the source galaxies, which we argue, could be furnished 
by polarimetric data. 



4.1 A spherical lens model 



To illustrate our argument we take an isothermal sphere with a core radius as our lens model (scheinder et al. 1992). The 
lens parameters are then given by 

k =!<&o(<r 3/2 + <r 1/2 ) (14) 

7 = i$ (d-3/2 „ d-i/2) (15) 

where d = (1 + (r/ro) 2 ) and r is the distance to the lens centre. In addition to the values of ro, and $0, the coordinates of the 
centre of symme try, introduce two further parameters. The centre of the lens can be inferred in various ways (see for example 
( tKochanek 1990|) ) and we shall assume in the rest of this article that it is already known. We have used a lens with ro = 1 



and $0 = 0.9 (Fig||). 



If we assume values for the pair of parameters (ro, $0), an d the image parameters are assumed to have been determined 
exactly, then it is possible to reconstruct the parameters for the source. If we assume also that the only information that 
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Figure 3. Lens parameters for the isothermal sphere, n long-dashed line, 7 short-dashed line, r = j/k full-line. 

we have about the source is from the measurement of the polarization, which provides us with the orientation of the source 
galaxies with a standard error of e, then we can find the lens parameters that best fit the observed polarization direction. 
The natural way to do this is to carry out a chi-squared minimisation, i.e. we minimise the following expression over the lens 
parameters. 



Xi(ro,$o) =2_^ -2 



(16) 



The summation in equation ( J16| ) is carried out over the set of galaxies whose polarization has been determined. a r a is the 
direction which has been reconstructed from the parameters of the lens model, and a™ is the measured direction. We have 
reconstructed lens parameters for the cases when the number of galaxies, N, for which the polarization has been measured is 
given by N = 50, 10 and 5. For each of these cases we consider the error on the polarization direction, e, to be given by 2, 5, 
and 10% .In order to get an idea of the confidence one should have in the values of ro and $0 obtained in this way, we have 
for, each pair of values (Nje), generated 500 different samples of galaxies, and from these used the chi-square estimator of the 
values of (ro,$o)- Figures ^| and ^| show ellipses in the (ro,"l , o) plane which contain 70% and 90% of the reconstructed values. 

In order to quantify the information that the measurement of the polarization provides, we have also carried out a 
reconstruction of the lens (that is determined ro and $0) assuming no knowledge of the orientation of the source galaxies. We 
did this by comparing the distribution of source parameters inferred from an assumed pair of lens parameters and sample of 
image galaxies with a distribution of source galaxies that is assumed to be uniform in orientation and ellipticity. We compared 
the joint distributions in ellipticity and orientation by chi-squared fitting of the two distributions. First we divide the ot s , e s 
plane (a s € [0, n] is the source orientation, and e s £ [0, 1] is its ellipticity ) into N c bins. We then minimise 



X2(ro,$o) 



Nc 

1=1 



n;(r ,$o) 



(17) 



with respect to (ro,$o)- n r is the observed number of galaxies in bin i and n % RX the predicted number for any assumed pair of 
lens parameters. We have binned the data in such a way that n l ex = 20. This procedure was carried out for several hundred 
samples of 500 galaxies. The reconstructed (ro,$o) are plotted on figure (^|). We can see that there is a great dispersion in 
the reconstructed data. This arises from the fact that the statistic of the reconstructed sources is fairly insensitive to the lens 
parameters. Therefore the above \ 2 is small for a wide range of parameters over which the minimum can be found. 

Although this method of finding confidence regions for the lens parameters in the absence of polarization information 
may not be optimal, it does however given some idea of the advantages of using polarimetric data (see Fig. 2) . 
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Figure 4. This figure shows the ellipses of the (rrj,$o) plane which contain 70% of the reconstructed points, the solid line is for 
[N = 50, e = 2%), the dotted for (N = 50, e = 5%), the small dashed for (N = 50, e = 10%), the long-dashed for (TV = 20, e = 2%), the 
dotted-small dashed for (N = 20, e = 5%) and the dotted-long dashed one for (/V = 10, e = 2%). 
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Figure 5. This figure is the same as Fig. [ij but with the ellipse containing 90% of the points. 
4.2 An arbitrary spherical lens 

As we have already mentioned, if one does not use the polarization to determine the orientation of the source galaxy, the only 
measurable quantities are the image parameters. We are thus left with five unknown parameters (Z„, Al 3 , a„, F, 9) but still 
have only 3 equations, and so it is impossible to solve the system. 

If we assume only that the lens is spherically symmetric, and, as before, that the centre of the lens is known, so that 
the angle 9 is determined, and that the polarization has been measured, then the system of equations (|lo|),(pT|) and ( |l3| ) will 
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Figure 6. Each point in this figure gives the couple (ro,&o) reconstructed from a sample of 500 galaxies without using the information 
given by the polarization. (There are no points outside this square because during the minimisation process, we limited the variation of 
r to [0.5, 2.1] and that of <J> to [0.4, 1.5] .) 
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Figure 7. The average error, expressed in percent, on the reconstructed value of T as a function of the error of the orientation, a s , of the 
source galaxy, expressed in degrees, One can see that if q s is given with a better than 5° precision then the error on T is less than 30%. 

contain three unknowns (l„, Al 3 , V). It is thus possible to deduce the remaining source parameters as well as F. From equations 
( [To[ ) and ( |l3| ) one can show that F is given by 



(ds a + SiC s )Ahr 2 + 2s s r + {as s - &c„)Afc = o . 



(18) 
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Figure 8. This figure shows the value and the error (square) on the parameter T as a function of radius. The bars are 95% error bars. The 
filled in triangles correspond to the values found for 10 randomly selected galaxies. The left (right) diagram corresponds to a population 
of galaxies with a 1 — a dispersion of 2.5° (5°) on the measure of a s . 

If there is no measurement error and if the lens is really spherical, this method would provide an exact determination of 
the lens parameters. Of course these assumptions are not very realistic, and it is interesting to investigate the case described 
by equation ([l8|) when there is a measurement error a s . In order to carry out this analysis, we have reconstructed T for 10000 
randomly selected galaxies as before, and including a random error on the measurement of a s . Fig.^ shows the mean error in 
percent on the value obtained for V (AT r — (T r — V) /Y where T r and Y are the inferred and the actual value) as a function 
of the error, measured in degrees, on a s (Aa s — — a s \ ). Evidently, if one can measure a 3 to within a few degrees, it is 
possible to locally determine I to a precision of around 20% using just one galaxy. 

We generated two samples of galaxies with a gaussian error distribution on a s with a standard deviation of 2.5° and 5°. 
The mean value and the one sigma error on the on the retrieved values V are shown in Fig |^. This figure shows that even 
with only ten or so galaxies for which the polarization hase been measured, it is possible to fairly accurately determine the 
function F(r). It is also important to point out that if one uses two galaxies at similar radial distances from the centre, it is 
possible to test the assumption of spherical symmetry. 

4.3 Arbitrary lenses 

In this section we discuss the problem of determining the lens parameters when no assumptions are made about the lens. If 
we no longer assume that the lens is spherical, then the angle 8 is not known a priori. The system of equations (pXj[) , ([Tl]) and 
( p^| ) then consists of 3 equations for 4 unknowns (l s , Al s , F, 8), and so it is impossible to solve for the unknowns. However if 
we assume that there are two galaxies in the sample that can effectively be associated with the same values of F and 8, then 
the previous system written for both galaxies will contain 6 unknowns, if, Al], Al^, T, 8), and six equations, which means 
that it can be solved exactly for T, 8 and the source galaxies' parameters. One should note, however, that if the orientation 
of the source galaxies is not known, adding more galaxies in this way cannot close the system, as there will always be more 
unknowns than independent equations: each new galaxy brings three unknowns and three equations, and so the equations 
are always underdetermined. There are a number of methods for finding galaxies associated with the same pairs of values of 
(T,8). If we make no assumptions about the lens, then one has to take two galaxies that have an angular separation much 
smaller than the angular diameter of the lensing cluster. One can then make the approximation that the potential is the same 
for both galaxies. If on the other hand we assume that the potential is spherically symmetric, then one can take galaxies 
which are well separated but at the same distance from the centre, so they have the same F and 9. There could of course be 
other symmetries that allow one to do this. 

The equation ( |li| ) for each galaxy gives the following system: 

/ (CfSJ 1 + S? 1 C«)Ai? i r 2 + 2STr + (Cf Si 1 - Sf CfOAZf =0 

\ {CpSI 2 + S? 2 C B /)AlfT 2 + 2S§ 2 T + (Cf Sf 2 - SfCf 2 )AZf = ( ' 

where g\ (resp. 52) denote quantities related to the first (resp. second) galaxy. The parameters Cf} and (k = i,s) depend 
only on 8. This system can be solved for F and 8. As before, we tested the possibility of reconstructing lens parameters 
by generating a sample of 10000 galaxies. Two galaxies are selected which correspond to lens parameters which are only 
slightly different. Thus we have considered only pairs of galaxies separated by less than 0.2 in the lens plane. To simplify the 
presentation we shall use the following notation. 
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Figure 9. Error on T plotted against the difference between T 91 and T 92 (see text) for galaxies for which ce s is measured exactly 
(squares) and a standard deviation of 2.5°. 

• r 9i parameter corresponding to the i* galaxy. 

• r retrieved value 

• T m = (r 91 + T 92 )/2, the average value. 

• Ari2 = |r 91 — r 92 |/r m : the initial difference in F for the two galaxies. 

• Ar — F — r m |/F m : the error on the retrieved value 

• Aq s = (Aq 91 + Aa 92 )/2. The quantities Aaf fc are defined as in the paragraph above. 

In Figj^ we have plotted the average value of Ar as a function of AFi2. If there is no measurement error, the error on 
the retrieved T is roughly equal to the difference between the values of T gi andr 92 . If on the other hand one plotsthe same 
curve but for a population of galaxies when a measurement error of 2.5° is made, then the error on V is noticeably greater, 
and around 15 — 20%. In order to establish the precision on a s necessary to establish F up to a certain level of accuracy, we 
have plotted the mean value of AF against Aa s for pairs of galaxies satisfying Ari2 < 10%. Indeed one sees that just one 
pair of nearby galaxies whose polarization yields a s with 5° accuracy allows T to be determined with an accuracy of 25% 
without making any assumptions about the form of the lens. Evidently, if one has a large number of galaxy pairs, that it is 
possible to still further reduce this error. 

One should also not that if one has several pairs of neighbouring galaxies, the number of equations in ( |l9| ) is greater than 
the number of unknowns, and one should use statistical techniques for the determination of the lens parameters. 



5 DETERMINATION OF k AND 7 

Once F has been determined, it is still necessary to determine the convergence ft and the shear 7, which are the parameters 
directly linked to the mass density distribution in the lens plane. These parameters can be obtained from equation (Q), which 
can be written in the form 



A 1 (i + r') + 2AA 1 rc* l 



(20) 



If we now average the above relation over a subset 5* of galaxies having the same value of k, either because they are 
sufficiently close together or because of some symmetry (all the galaxies can be used and not only those with a measured 
polarization)we obtain k 2 = As/a, where a, the average of a over the subset S, can be computed using the image and the 
reconstructed value of F and 9 (which are equal for all galaxies in S). The average of A s , A s , cannot be computed since the 
intrinsic brightness of the sources are not known, but since the source galaxies (and therefore their properties) are uniformly 
distributed in the source plane, A s should be constant (independent of the set 5* chosen) if the average is taken over a 
sufficiently large number of galaxies. Therefore we can arbitrarily set A s = l for the whole subset of galaxies. This allow to 
compute ii = l/o and 7 = Tk which are proportional to k and 7 respectively, the proportionality factor being a constant 
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Figure 10. Error on T plotted against the average error in a 3 . 

throughout the lens plane. The reliability of this method will strongly depend on the number of galaxy contained in each 
subset and on the luminosity dispersion of the source galaxies. Each subset should contain enough galaxies in order for to 
assume the A s is equal to the average luminosity of the source galaxies. 



6 OBSERVATIONAL CONSIDERATIONS 

In the previous section we saw how knowledge of the orientation of the source galaxy enabled us to improve the determination 
of the gravitational field of the cluster. In principle this orientation would be given by the direction of the polarization. 
Similarly measurement of the degree of polarization would yield information about the ellipticity of the source, although in 
view of the possible uncertainty concerning the mechanisms involved in the production of the polarization detailed modelling 
might be required. Potentially, determination of the intrinsic polarization of lensed galaxies could yield valuable information. 
Thus a question of crucial importance is the feasibility of such measurements and the reliability of their interpretation. 

The main difficulty evidently arises from the faintness of the source galaxy, and the relatively low level of polarization. 
It should be borne in mind, however, that although we are looking at objects of magnitude 25 or fainter, it is the integrated 
polarization that we are interested in, and measurements could be taken over broad band. Such observations in the optical 
have been carried out for a number of high redshift radio galaxies (Tadhunter et al. 1992), with the galaxies showing high 
polarization, which appears to increase with higher redshift. Although our concern is mainly with spiral galaxies, it is possible 
that radio galaxies could also be used in the analysis, since these appear to display net polarization oriented in a definite 
direction to the radio axis. Our main point here however is that optical polarimetric measurements would be difficult, but not 
impossible, particularly on a large telescope. It is important to point out also that in our analysis it is only the direction of 
polarization that is important, and this is much easier to establish accurately than the degree of polarization (Treanor 1968). 

Another problem would be the possible effect of Faraday rotation, yet in the visible we would expect this to be fairly 
weak. If we were to use polarimetric measurements in the radio, then this of course would need to be taken into account. 



7 CONCLUSIONS 

The polarimetric study of weakly lensed galaxies can in principle provide information about the lensing gravitational potential 
that otherwise would be extremely difficult or impossible to obtain. The method relies on the property that the plane of 
polarization of light is not affected by the gravitational potential. Thus the direction of polarization provides a marker for the 
orientation of the source galaxy, and hence information about the distortion in the image induced by the lens. This reduces 
the uncertainty in the estimation of the lens parameters, and can even, in certain situations when a few source galaxies can 
actually be measured polarimetrically, break the degeneracy of the equations that determine the lens parameters. 

Although to a large extent we have assumed that Thomson or Rayleigh scattering is responsible for the polarization flux, 
the our conclusions are not greatly affected by the presence of such small scale mechanisms as dichroism provided that on 
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the large scale axial symmetry holds. Only in exceptional circumstances would one expect the direction of the polarization 
not to be perpendicular to the major axis of the elliptical isophotes. Since weak lensing would induce a rotation of much less 
than 90 degrees, even if the integrated polarization were aligned parallel to the major axis of the source, this should not lead 
to confusion and hence measurement of the direction of polarization would still provide rigorous constraints on the lensing 
parameters. 

Measurement of the optical polarization of faint galaxies would require large telescopes and integration times, but it is 
important to realise that the main advantage of the proposed method relies on the use of the direction of polarization, and 
not the degree of polarization and the former is considerably easier to establish. 
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